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The transition of a system of partial differential equations which describe the sta- 
tionary flow behind the shock-wave front of a detonation complex upon detonation 
of a cylindrical charge to a system of ordinary differential equations is performed by 
means of the series expansion in terms of the radial variable. The necessary equations 
for determination of the derivatives of solutions with respect to the parameters and 
the initial conditions for them are formulated. Imposing the condition of continuous 
extendibility of the solutions leads to equations that allow one to determine the shape 
of a shock-wave front as a function of wave velocity. 



The attempt is made to construct a consistent 
method of determining the shape of the shock-wave 
front of a detonation complex upon detonation of a 
cylindrical charge of finite diameter in the case where 
the detonation velocity deviates insignificantly from 
the Chapman- Jouguet velocity. This problem can 
be solved, disregarding the problems connected with 
the description of the flow near the charge edge. This 
problem was considered earlier (see, e.g., [1, 2]), but 
in the solution of it, the authors restricted them- 
selves only to the determination of the curvature of 
the front, i.e., the first expansion coefficient of the 
function describing the front shape in terms of the 
powers of its argument. The distinguishing feature 
of the proposed approach is the fundamental possi- 
bility of determining the shape of the front with no 
matter how high accuracy. 

We write a system of equations that determines 
the axisymmetric stationary flow in the presence of 
a chemical reaction characterized by the unique vari- 
able A (fraction of the reacted substance) in the cylin- 
drical coordinate system: 
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Here u and v are the radial and axial velocity-vector 
components, respectively, E is the internal energy per 
unit mass, c is the velocity of sound with constant A, 
and R is the chemical-reaction rate which depends, 
in the general case, on the pressure p, the density p , 
and the variable A (for brevity, the arguments of this 
function are omitted). 

We pass to the coordinate system (l,r) deter- 
mined by the relation l = z — Zf(r), where Zf(r) is 
a function describing the front surface in cylindri- 
cal coordinates. In coordinates (Z,r), the curvilinear 
front is a plane 1 = 0, which simplifies the writing of 
the initial boundary conditions. In these coordinates, 
system (1) is reduced to the form 
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{ V + ud)— + u — = R , 



where d = —z'f(r). 

With the aim at finding the shape of the front at 
a given velocity of the detonation wave, we expand 
the function Zf(r) in powers of r: 



OO 

z f( r ) = Z a * r2j; 



other functions: 

OO 

P(r,l) = ^2pi{l)r 2 \ 

i = 0 



P(r,l) = ^ 2 pi(l)r 2 \ 
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v(r, l) = ^2 v i(l) r2 \ 

i = 0 



A(r,0 = ^A,(/)r 2i , 
2—0 



u(r,l) = y^Uj(/)r 22+1 . 
2—0 



[Only the even powers of r in the expansion of p, p , 
v , and A and only the odd powers in the expansion of 
u follows from the requirement for the analyticity of 
these functions on the (r = 0) axis. We also note that 
the functions with zero subscripts are the values of 
the corresponding quantities on the axis, except for 
uq{1) equal to du/dr for r = 0]. Substituting these 
series into system (2) and grouping the coefficients 
at the same powers of r, we obtain the system of 
ordinary differential equations for the functions pi{l), 
Pi(l ), Vi(l), Ui(l), and A* (Z): 
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Here the prime denotes the derivative with respect to 
the variable l; R m (l), K m (l ), and c^(7) are the series 
expansion coefficients of the following functions in 
terms of the powers of r: 
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We consider the dependence of the solutions of 
this system on the coefficient a k and the detonation- 
wave velocity I? as the parameters. Resolving the 
equations of system (3) relative to the derivatives, 
one can note that, generally, the solutions of these 
equations can be extended only up to the point 
i>o = cq- We write, for example, the normalized 
(i.e., resolved relative to the derivatives) equations 
for to = 0: 
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(The equations are similar for other values of to.) 
The solutions of these equations are extendible 
through the point Vq = cq only if the expression 
2cqPqUo — K 0 vanishes simultaneously with the ex- 
pression — eg. Precisely the condition of continu- 
ous extendibility of solutions through this point (i.e., 
the absence of a singularity) determines the values of 
a k (shape of the front) as a function of D. 

To find the dependence a k (D ), we do the follow- 
ing. Let the coefficients a k be known at a certain 
velocity D (a k = 0 for D = -Dcj, where Dqj is the 
velocity of the Chapman- Jouguet ideal detonation 
in an unbounded medium, i.e., in the case of a plane 
front). We find the derivatives da k /dD , d 2 a k /dD 2 , 
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etc., thus determining the dependence a k (D). To do 
this, we note that the derivative 

dpm _ d'Prn dp m da k , ... 

dD ~ dD + ^ da k dD 1 j 

k 

(here one can use also the density or velocity deriva- 
tive instead of pressure derivatives) is finite at the 
point Mq = Uo/co = 1, whereas the derivatives 
dpm/dD and dp m /da k tend to infinity in approach- 
ing this point; this can be seen after normalization 
of the differential equations for these functions (see 
below). Equality (4) also allows us to find da k /dD. 
We find these derivatives for a k = 0 and D = Dcj. 
We write 

dp 0 _ dpo d]Po_ da i dpo da 2 
dD ~ dD + da x dD + da 2 dD + ' ' ' ’ 



dpi dpi dpi da x dpi da 2 

dD dD da x dD da 2 dD ~*" 

etc. As is shown below, dp m /da k and other deriva- 
tives are different from zero only for in < k, and 
dpm/dD is different only for m = 0. Whence, con- 
sidering the limit Mq — > 1, we find that da k /dD = 0 
for k > 1 and 

da l = _ lim dPo/dD (5] 

dD Mo — >1 d Po /d ai ’ 1 ’ 

i.e., in the D-linear approximation, only the first coef- 
ficient of the expansion Zf(r ) in terms of the powers of 
r (curvature of the front) depends on the detonation- 
wave velocity. We consider the next approximation 
to find d 2 a k /dD 2 . By analogy with the previous case, 
we write the derivatives d 2 p m / dD 2 , discarding, with 
allowance for the aforesaid, the zero terms: 

d 2 po _ d 2 p 0 d 2 p 0 (dai\ 2 0 d 2 p 0 dai 
dD 2 “ dD 2 + daj V dD) + da x dD dD 

dpo_ d 2 a x dpo d 2 a 2 dpo_ d 2 a 3 
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d 2 p 1 _ d 2 pi d 2 pi / da\ \ 2 0 d 2 p x da x 

dD 2 ~ dD 2 + da 2 V dD ) + da x dD dD 
dpi d 2 a 2 dpi d 2 a 3 

da 2 dD 2 da 3 dD 2 

d 2 p 2 _ d 2 p 2 d 2 p 2 / da x \ 2 0 d 2 p 2 dai 

dD 2 ~ dD 2 + da 2 V dD ) + da 3 dD dD 

dp2 d 2 a 3 
da 3 dD 2 

etc. We show below that the derivatives d 2 p m /da\ 
and d 2 p m /daidD are zero for to > 1, and d 2 p m /dD 2 
is zero for to > 0; therefore, considering again the 
limit Mo — > 1, we obtain d 2 a k /dD 2 = 0 for k > 2 
and 
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( 7 ) 



Thus, to calculate da x /dD, it suffices to know 
how to find the derivatives of the functions p m with 
respect to the parameters a k and D (in this case, the 
functions dpo/da x and dpo/dD), and, in addition, 
the functions d 2 po/dD 2 , d 2 po/da\, d 2 po/da\dD 1 
dpo/da 2 , dpi/da 2 , d 2 pi/da\, and d 2 pi/da x dD to 
calculate d 2 ai/dD 2 and d 2 ai/dD 2 . 

Similarly, one can, in principle, find any deriva- 
tive of the coefficients a k with respect to the detona- 
tion velocity D, i.e., to determine the shape of the 
front with required accuracy. 

We now consider the procedure of calculation of 
the above-mentioned derivatives of the functions p m 
with respect to the parameters a k and D entering 
formulas (5)-(7). To start with, we note that we 
are interested in the derivatives at a fixed value of 
Mo rather than at a fixed value of l; equations for 
the latter can be derived by varying the equations of 
system (3) relative to the parameters a k or D (see, 
e.g., [3]). The relationship between these derivatives 
has the following form: 

dpm{M 0 ,a,D) = dp m (l, a, D) 
da k da k 

dpm(l,a,D) dl(M 0 ,a,D) 
dl da k 



dl(M 0 ,a,D) dMo(l,a, D) /dM 0 (l,a,D) 

da k da k / dl 

Here a denotes the set of coefficients a k , where k = 
1, . . . , 00 . A similar relationship between the second- 
order derivatives has the form 
d 2 pm{M 0 ,a,D ) _ d 2 p m (l, a, D) 
daidaj da x daj 

d 2 p m (l, a, D) dl(M 0 ,a,D) 
daidl daj 

d 2 p m (l, a, D) d?(Af 0 ,a,T>) 
dajdl dai 

d 2 p m (l,a,D) dl(Mo,a, D) dl(Mo,a,D) 

+ dH do) d^j 

dp k (l, a, D) d 2 l(M 0l a 1 D) 
dl daidaj 
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d 2 l(M 0 ,a, D) sd 2 Mo(l,a,D) 

dctidcij V dciidcij 

d 2 Mo(l,a, D) dl(M 0 , a, D) 
ddidl ddj 

d 2 M 0 (l,a, D) dl(M 0 ,a, D) 
ddjdl ddi 

d 2 Mo(l,a,D) dl(Mo,a,D) dl(M 0 ,a,D)\ 

+ m ha, d^i ) 

/dM 0 (l,a,D) 

/ dl 

As has already been mentioned, one can find the 
equations for determination of the derivatives of the 
functions p m (l, a, D) with respect to the parameters 
dk or D by varying the equations of system (3) rel- 
ative to the parameters a k or D. We write, for ex- 
ample, the system of equations for determination of 
dpm/ddk for a k = 0 and D = D c j [here p m (l) = 0, 
p m (l) = 0, etc., for to ^ 1 and u m {l) = 0, and the 
functions with zero subscripts are corresponding ex- 
act functions in the plane case, i.e., upon detonation 
in an unbounded medium]: 
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Here 

d K m _ dK_ d Pm dK_ dp m dK_ dX m 

ddk dp ddk dp ddk d\ ddk ’ 
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The derivatives with respect to p, p, and A are taken 
for p = po(l), P = Po(P), and A = Xo(l), respectively. 
Considering this system with the initial conditions 
(see Appendix), one can conclude that dp m /ddk and 
the other derivatives are different from zero only for 



to < k (by virtue of the fact that m ^ k, the ini- 
tial conditions are zero, and the equations are homo- 
geneous). The equations for dp m /dD can be writ- 
ten in the same way; their consideration allows us to 
draw a conclusion that <9p m /ctD = 0 for to ^ 0, and 
du m /dD = 0 for any value of to. 

We also present the system that determines 
d 2 Pm /ddf f° r a k = 0 and D = Dcj with allowance 
for the fact that only the derivatives dpo/ddi, 
dpo/ddi, dvo/ddi, duo/ddi, and dXo/dd\ are not 
zero: 
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Here 
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2 , d2R °> 
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dp 2 1 
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dp 2 ' 
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d 2 R 0 /dA 0 \ 2 d 2 R 0 dp 0 dp 0 
dX 2 V da\ ) dpdp da\ dci± 
d 2 R 0 dp 0 dX 0 d 2 R 0 d\ 0 dp 0 

dpd X da\ da\ d\ dp da\ da\ 

dR d 2 po dR d 2 po dR d 2 Xp 
+ ~d p ~fa 2 + ~dp ~fa[ + ~dX "da 2 ’ 

d 2 R m _ dR d 2 p m dR d 2 p m dR d 2 X m 
da 2 ~ dp <9a 2 + <9p da 2 + dX da 2 

for to > 0, 

d 2 K m /da\ and <9 2 c^/<9a 2 have a similar form. One 
can conclude from the consideration of this system 
with the boundary conditions that d 2 p m /da\ = 0 
for to. > 1. Similarly, we have d 2 p m / da\dD = 0 for 
to > 1 and d 2 p m /dD 2 = 0 for to. > 0. 

Equations (8) and (9) and the initial conditions 
formulated for them in the Appendix allow us to find 
the derivatives of the solutions with respect to the 
parameters a k and D , which, in turn, allow us to 
find the shape of the shock-wave front in the form 
of an expansion of the coefficients a k in terms of the 
velocity D. A similar procedure can be performed to 
find the derivative a k with respect to D of any order; 
however, it is necessary to note that the complex- 
ity of calculations rapidly increases with the order of 
derivatives. 

Thus, the problem of determining the shape of a 
detonation wave front as a function of the detonation 
wave velocity upon detonation of a charge of finite di- 
ameter has been solved. In addition to the fact that 
this problem is interest in itself, it is a necessary com- 
ponent of the solution of the more general problem of 
finding a relation between the basic quantities which 
characterize the detonation of a cylindrical charge: 
the detonation wave velocity, the radius of a charge, 
and the shape of a wavefront. 

The author thanks N. M. Kuznetsov and V. G. 
Grudnitskii for valuable comments and discussions. 



APPENDIX 

We present the values of the density, pressure, 
and velocity components behind the front of a curved 
shock wave in the form 

Pf(r, D, a) = p id (D cos cc), 

Pf(r,D,a) = Pi d (D cos a), 

Vf(r,D,a ) = Dcos^cosabtan 2 a + p 2 / p 2 d , 



Uf(r,D,a ) = I? sin <pcos cc^/tan 2 a + p 2 / p 2 d , 

X (r, D, a) = 0, 

(p = arctan (tan a(pid/pb)) — a, 

where a is the angle between the normal to the front 
at a given point and the direction of detonation-wave 

OO 

propagation, tan a = z'Ar ) = ^ 2ia i r 21 ^ 1 , and p id 

i = 1 

and p^ are, respectively, the density and the pressure 
behind the front of a plane shock wave propagating 
with the velocity D cos a in the same initial medium, 
i.e., at the same values of pb and pb (the density and 
the pressure in front of the shock wave). The initial 
conditions for system (3) have the form 
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r—0 



r=0 



(2to + l)!<9r 2m+1 1 

The initial conditions for system (8) have the follow- 
ing form: 
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We modify their form. For example, for the pressure 
and the radial velocity component, we write 



&Pm 
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dp f 


d tan a 


da k 


i=o (2 m)\dr 2m 


d tan a 


da k 



ai—0,r=0 



du„ 

da k 



^2m+l 

i=o (2m + l)\dr 2m+1 



duf 


d tan a 


d tan a 


da k 



ai=0,r—0 



Since pj is an even function of tan a, and a = 0 
for a* = 0, we have (dp m /da k )\i = o = 0. For 
the same reasons, we have (dp m /da k )\i=o = 0 and 
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( dv m /dak)\i=o = 0. The velocity Uf is an odd func- 
tion of a; therefore, 



d u f 


(9 tan a 


d tan a 


da k 



i.e., only (duk-i/da k )\i=o are different from zero. 
Similarly, since 

d 2 p d tan a d tan a 
d tan a 2 dat dcij 



r 2i— l r 2j — 1 



only (dp i+ j-i/daidaj)\i= 0 , (dp i+j - 1 /da l da j )\i =0 , 
and (dvi+j-i/daidaj\i—o are different from zero. 
Using similar considerations, one can show that 
among the derivatives d/dD , only the derivatives 
dpo/dD, dpo/dD, dv 0 /d D, d 2 p 0 jdD 2 , d 2 p 0 /dD 2 , 
d 2 Vo/dD 2 1 and d 2 uo/dDdai are not zero. 
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